Transport Gap in Suspended Bilayer Graphene at Zero Magnetic Field 
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We report a change of three orders of magnitudes in the resistance of a suspended bilayer graphene 
flake which varies from a few kQs in the high carrier density regime to several MQs around the 
charge neutrality point (CNP). The corresponding transport gap is 8 meV at 0.3 K. The sequence 
of appearing quantum Hall plateaus at filling factor u = 2 followed by = 1 suggests that the 
observed gap is caused by the symmetry breaking of the lowest Landau level. Investigation of the 
gap in a tilted magnetic field indicates that the resistance at the CNP shows a weak linear decrease 
for increasing total magnetic field. Those observations are in agreement with a spontaneous valley 
splitting at zero magnetic field followed by splitting of the spins originating from different valleys 
with increasing magnetic field. Both, the transport gap and B field response point toward spin 
polarized layer antiferromagnetic state as a ground state in the bilayer graphene sample. The 
observed non-trivial dependence of the gap value on the normal component of B suggests possible 
exchange mechanisms in the system. 
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I. INTRODUCTION 

Followed by the isolation of single layer graphene, the 
study of bilayer graphene (BLG) became a separate di- 
rection of research in the community of two dimensional 
materials. Charge carriers in bilayer graphene have a 
parabolic dispersion with an effective mass of about 
0.054mefii^ but also possess a chirality. The latter man- 
ifests itself in an unconventional quantum Hall effect^ 
with the lowest Landau level (LLL) being eight fold de- 
generate. Compared to single layer, bilayer graphene 
has next, to spin and valleys degrees of freedom, an ad- 
ditional orbital degree of freedom, where Landau lev- 
els with numbers n = and 1 (each four fold degen- 
erate) have the same energyi^ Recent advances in ob- 
taining suspended bilayer graphene devices with charge 
carrier mobility exceeding fi > 10,000 cm^V~^s~^ gave 
access to the investigation of many-body phenomena in 
clean bilayer graphene at low charge carrier concentra- 
tion (n < 10^° cm-2)4-ii 

Due to the non vanishing density of states at the charge 
neutrality point (CNP), bilayer graphene is predicted to 
have a variety of ground states triggered by electron- 
electron interaction. There are two competing theories 
describing the ground state of BLG: a transition (i) to 
a gapped layer polarized state (excitonic instability 
or (ii) to a gapless nematic phaseJ^"— 

Excitonic instability is a layer polarization in which 
the charge density contribution from each valley and 
spin spontaneously shifts to one of the two graphene 
layersji^iii This redistribution is caused by an arbitrar- 
ily weak interaction between charge from conduction and 
valence band statesJ^ii^ Since each bilayer flavor (spin 
or valley) can polarize towards either of the two layers. 



there are 16 possible statesfi^iii which can be classified 
by the total polarization as being layer ferromagnetic (all 
degrees of freedom choose the same layer), layer ferri- 
magnetic (three of the four valley-spin flavors choose the 
same layer), or layer antiferromagnetic (with no overall 
polarization). To make it clear, the therm "magnetic" 
should be associated to flavors (not only spin) orienta- 
tion in between two layers. These states are considered 
as analogous to the biased bilayer— in the sense that 
the charge transfer can be attributed to the (wave vector 
dependent) exchange potential difference between low- 
energy sites on the opposite layers The total energy of 
the system is lowered by the gain in the exchange inter- 
action via breaking of the inversion symmetry, i.e. intro- 
ducing a gapped state. Antiferromagnetic polarization 
is electrostatically favorable due to the absence of a net 
charge on both layers, however, the actual ground state is 
theoretically undefined J^ii^iii Recent experiments have 
suggested the evidence of the possible existence of two 
of the antiferromagnetic states - the anomalous quan- 
tum Hall state (AQH)^i^ and spin polarized layer anti- 
fcrromag netic state (LAF)1. To avoid possible confusion 
we note that in earlier literatures^ the LAF state is also 
called quantum valley Hall state. The AQH has electrons 
being polarized in the same layer for both spins and in 
opposite layers for opposite valleysJ^iii This state has 
spontaneously broken time reversal symmetry and there- 
fore possess a substantial orbital magnetization exhibit- 
ing quantized Hall effect (at zero magnetic field), while 
its spin density is everywhere zero Due to its magneti- 
zation the AQH can be favored over other ground states 
in the perpendicular magnetic field. The LAF state has 
opposite spin-polarization for opposite layers. In con- 
trast to AQH, the LAF state does not have topologically 
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protected edge states, which brings its minimum conduc- 
tance to zero. For both states the theoretical estimations 
of the gap A give the value of 1.5-30 meVM^ How- 
ever, the inter-valley exchange weakly favors the LAF 
stateJ^i^ One of the ways to determine the character 
of the bilayer ground state experimentally is to investi- 
gate the response of the gap value to the magnetic field 
B (which couples to spin) and electrical field E (which 
couples to layer pseudospin) ^ When Zeeman coupling is 
included, the QAH state quasiparticles simply spin-split, 
leaving the ground state unchanged but the charge gap 
reduced. It was calculated that for a 4 meV spontaneous 
gap at zero-field, a field of B = 35 T drives the gap to 
zero. On the other hand, the gap of LAF is weakly B 
field dependent. 



The second possible description for the ground state 
of BLG is based on a nematic phase caused by the 
renormalization of the low energy spectrum i^^i^^ Detailed 
tight-binding model studies showed that including next- 
neighbor interlayer coupling changes the band structure 
in bilayer producing a Lifshitz transition in which the 
isoenergetic line about each valley is broken into four 
pockets with linear dispersioni^i^ At the energies higher 
then 1 meV the four pockets merge into one pocket 
with the usual quadratic dispersion. Moreover, electron- 
electron interactions might result in the further energy 
spectrum transformation, where the number of low en- 
ergy cones can be reduced to two near each of the two K 
points.— In this case the minimum conductance of the 
bilayer graphene is supposed to be increased comparing 
to bilayer with parabolic dispersion (8e^//i). This sce- 
nario was also supported by the experimental result on 
the suspended bilayer graphene in which strong spectrum 
reconstructions and electron topological transitions were 
observed 3 



In this paper we present electric transport properties 
of suspended bilayer graphene by studying its behavior in 
tilted magnetic fields. At B = T we observe the sponta- 
neous opening of a gap by changing charge carrier density 
from the metallic regime (n = 3.5 x 10"'^"'^ cm~^) to the 
CNP. At a temperature of 1 K we measure a resistance in- 
crease from 5 kfi up to 14 Mfi. The observation indicates 
the gapped ground state of the studied bilayer graphene 
with a value of 6.8 meV. Measurements in tilted magnetic 
field showed that the resistance at the CNP decreases 
with an increasing of magnetic field. Based on this we 
propose a possible scenario of the symmetry breaking in 
this bilayer graphene sample: Spontaneous valley split- 
ting at zero magnetic field followed by the splitting of 
the spins originating from different valleys with increas- 
ing of B. Both, the gap value and its weak linear decrease 
with B, supports LAF as the ground state of the studied 
sample. 



II. EXPERIMENTAL DETAILS 

Suspended bilayer graphene devices were prepared us- 
ing an acid free technique.— We deposited highly or- 
dered pyrolytic graphite on n"'""''Si/Si02 wafer (500 nm 
thick) which is covered with an organic resist LOR 
(1.15 /xm). A standard lithography procedure is per- 
formed in order to contact bilayer graphene flakes (de- 
termined by their contrast in optical microscope) with 
80 nm of Ti/Au contacts. A second electron beam lithog- 
raphy step is used to expose trenches over which graphene 
membrane becomes suspended. To achieve high quality 
devices we use current annealing technique by sending a 
DC current through the membrane (up to 1.1 mA) at the 
temperature of 4.2 K. While ramping up the DC current, 
simultaneously, we keep track of the sample resistance. 
Once the resistance reaches values in the order of 10 kfis 
we stop annealing and check the gate voltage dependence. 
We repeat this procedure till the appearance of a sharp 
resistance maximum at the CNP located close to zero 
Vg. More details on the current annealing procedure 
can be found in Tombros et ali^ The studied device was 
2 /im long and 2.3 /im wide. All measurements were per- 
formed in four-probe geometry (with contacts across the 
full width of graphene) at the temperatures from 4.2 K 
down to 300 mK. The four-probe method allows to elim- 
inate contact resistances. As discussed below the resis- 
tance measurements consist of a superposition of longitu- 
dinal magnetoresistancc {pxx) and Hall-resistance {pxy)- 
The carrier density in graphene is varied by applying a 
DC voltage (Vg) between the back gate electrode [n^'^ Si) 
and the graphene flake. Based on serial capacitors model 
a unit capacitance of the system is 7.2 aF^um^^, which 
relates gate voltage with density as n = aVg where a is 
a leverage factor of a = 0.5 x 10^°cm~^V~^. The typical 
current we use is around 1 nA. See Appendix. 



III. TEMPERATURE DEPENDENCE AND 
QUANTUM TRANSPORT 

Our pristine samples are strongly p-doped with the 
CNP situated beyond 60 V and a metallic resistance of a 
few hundreds of Os over the entire voltage range. There- 
fore we perform current annealing technique in order to 
obtain high quality devices. In contrast to previous sam- 
ples, in which each next step of current annealing tend 
to cause sharper change in the resistance values within 
the scanned region of Vg, the discussed bilayer sample 
already shows after the flrst current annealing step a 
high resistive region around the CNP (not shown). The 
next annealing step (1.1 mA) moves the charge neutrality 
point down io Vg = i V. However, surprisingly the re- 
sistance around CNP becomes 14 MJl and reduces down 
to 5 kfi in the metallic regime at Vg = -60 V (Fig. [Iji, 
Inset). This fact points toward opening of a gap. The 
temperature dependence of the membrane from 4.2 K to 
300 mK is shown in Fig. [T^. There is an essential change 



3 




-15 15 
V (V) 



0.4 0.6 
1/T(1/K) 



FIG. 1. (Color online) 4-Probe resistance of the suspended 
bilayer graphene. a) Gate dependence of the sample at the 
temperatures of 4.2 K (in black), 1.3 K (in red), 0.7 K (in 
green), 0.3 K (in blue). Inset: Resistance at 4.2 K in log- 
scale showing the dramatic change from the CNP to metallic 
regime; b) Transport gap extraction at 4.2 and 0.3 K. The 
energy gap in bias direction is highlighted by the conductance 
crossover (fitted with dashed lines) at zero. The values of the 
transport gap are 3 meV (4.2 K) and 8 meV (0.3 K); c) An 
Arrhenius plot of the resistance. The value of an extracted 
thermal gap is 0.33 meV 



of about 6 MJ7 in the maximum resistance {Rmax) from 
4.2 K down to 1.3 K, however further lowering of temper- 
ature does not change Rmax much. From an Arrhenius 
plot of the resistance at CNP (Fig. [T]:) we can extract 
a thermal excitation gap of 0.33 meVi^ The saturation 
of resistance at lower T can be explained by a variable 
range hopping with different temperature dependence. 
We would like to point out that our excitation current 
value of 1 nA gave a voltage drop of cx 10 mV at the 
CNP, which is much higher than kT energy at measured 



temperatures (0.3 meV). Therefore one has to be careful 
in comparing transport and thermal excitation gaps. 

There might be a couple of scenarios for the observed 
gap formation in the gate voltage dependence: (i) A lat- 
eral confinement in membrane, where energy levels are 
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W = 2.2) fiTti - width of the flake, I is integer value. How- 
ever, first two levels have energies oi Ei ~ 1.3 /.jcV and 
E2 ~ 5.3 /leV, which is much lower than /cgT at measured 
temperatures, (ii) True gap formation with zero density 
of states within the gap and available states at the con- 
duction and valence bands, (iii) Transport gap, accompa- 
nied by the observation of the reproducible conductance 
oscillations in the region of suppressed conductance. In 
such regime transport is limited by the quantum confine- 
ment effect along the width (mainly originating from the 
impurities) jiSI (iv) More complicated case, when the gap 
value depends on the charge carrier density, i.e. the en- 
ergy of the levels changes while being filled with carriers. 
This situation might happen when the gap is induced 
by charge redistribution in between layers, which would 
be influenced by the applied back gate voltage. At the 
moment, we can not determine the exact gap type, there- 
fore, further analysis is performed assuming a transport 
gap scenario, but keeping in mind that this gap value can 
depend on the density. 

In an analogy to graphene nanoribbon studies^li^^ we 
extract the transport gap from the gate dependence of 
the sample conductance as shown in Fig. [TJd). From a 
linear approximation of conductance one gets a region 
of AVg, where sample shows insulating behavior. This 
region AVg relates to the wave vector as Afc = VTrAn = 

TraAVg. Taking into account the quadratic dispersion 
of bilayer graphene, the corresponding energy scale can 
be calculated as 



AEf = —— = —naAVg 
2m 2m 



(2) 



From conductance graphs at different T we find 
AEf = 3 meV at 4.2 K and AEp = 8 meV at 0.3 K. 
The values of the transport gap are comparable to the 
energy gap (extracted in bias direction) values of sin- 
gle layer graphene nanoribbons of 50-85 nm wide^^i^ 
where in contrast to our case the gap is created by lat- 
eral confinement. The resistance value of 5 kfJ in metallic 
regime, similar to regular graphene devices, serves as an 
additional justification of excluding a lateral confinement 
as a cause of the observed transport gap. We can calcu- 
late the mobility of the charge carriers using a standard 
formula fi = l/(ei?sqn), where Rsq is a square resistance 
of the sample and e is elementary charge. The mobility 
value fj.(x 20, 000 cm^ V"^s"^ at n = 3.5 x lO" cm'^ cor- 
responds to the value of high quality bilayer graphene de- 
vices. Due to the symmetry of resistance change around 
CNP (Fig. [TJd) and the fact that CNP itself is situated 
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Given the fact that the resistance values reach Mfis at 
the CNP, it is already hard to establish quantum Hall 
plateaus in our suspended bilayer device. However, we 
have achieved to observe quantum Hall transport shown 
in Fig. [2ji,b). First quantum Hall plateau appears at 
5 T on electron side (red curve), which we attribute to 
the filling factor v ^ 2. This plateau is followed by the 
appearance of = 1 at 7 T (Fig. [^b)- The conductance 
values of the observed plateaus deviate from the expected 
ones of 2e^ /h and le^/h, since they are affected by charge 
inhomogeneity. Therefore, we determine the exact values 
of the corresponding plateaus by the scaling of their po- 
sitions in density {Vgi,) with magnetic field B (Fig. [2jl). 
As expected from i/ — n/{eB/h) the scaling is linear with 
the leverage factor of a = 0.64 x 10^°cm~^V~^ for i/ = 2 
and 1. In order to use the same a for both filling fac- 
tor sets (see Fig. [2}1) the slopes of versus B; and i/ 
values respectively, have to be twice as different. There- 
fore, we have to point out that the linear scaling will 
hold as well for a leverage factor of 1.1 x lO^^cm^^V^^ 
in case we assume v = A and 2 as an observed sequence of 
plateaus. From previous studiea^S we know that capaci- 
tance probed by the QHE in graphene devices (especially 
in suspended samples) can be higher than the geometri- 
cal value, due to the deviation from the plane capaci- 
tor model. However, we attribute the observed plateaus 
to the filling factors 2 and 1. As we noticed before^^ii^ 
most of the time the current annealing procedure leads 
to the formation of high quality annealed regions con- 
nected in series with low mobility p-dopcd regions close 
to the contacts. Therefore higher values of the conduc- 
tance plateaus can be explained by a "p-doped" slope, 
which increases with magnetic field B. This can be also 
the reason of the absence of resistance quantization in the 
electron-side (Fig. [2t). Assuming fiB ^1 for the forma- 
tion of QHE plateaus,— our observation implies a lower 
bound for the mobility of 2,000 cm^V~^s~^. 



FIG. 2. (Color online) Quantum transport at 1.3 K. a) Quan- 
tum Hall conductance of the suspended bilayer at zero and 
B = 5 T; b) Quantum Hall conductance of the suspended 
bilayer at B = 7 and 11 T. The exact filling factors v corre- 
sponding to the observed plateaus are shown; c) Resistance 
of the sample in quantum Hall regime, d) Scaling of the fill- 
ing factors position in gate voltage {Vgv) with magnetic field; 
e) LL hierarchy of the symmetry-breaking of the lowest LL 
in bilayer graphene. Suggested scenario of the spontaneous 
valley splitting followed by the spin splitting at high B. 



around zero gate voltage {Vg ~ 1.2 V), that corresponds 
to the density of n = 0.77 x 10^" cm~^ at V, we can 
also exclude the low quality "p-doped" regions close by 
the contacts (which can form after current annealing) as 
the cause of the reported gap. In the meantime, we can 
not exclude a charge inhomogeneity in the sample bulk 
which might lead to the observed order of magnitude def- 
erence between electrical and transport gaps, in analogy 
to nanoribbon case. 



To summarize our QH transport results: At this point 
we have shown that a zero-field gap opens at the CNP in 
the studied graphene bilayer. This observation points out 
on a possible symmetry breaking of the ground state in 
bilayer graphene. The application of B does not restore 
the broken symmetry and brings the systems in to the QH 
regime. In Fig.[2t) we show the hierarchy of the splitting 
of the 8-fold degenerate lowest Landau level in applied 
B^ The development of the level structure with B will 
be specified and discussed in section 4. In the meanwhile, 
if we assume that at _B = T one of the degeneracies is 
already lifted, then, with increasing field, one can expect 
quantization at the filling factors i> = and 4 followed by 
v = 2 and 1. However, if the initial symmetry breaking 
is strong enough and the scanned window in energy is 
limited (Vg), then one can expect quantization at v = 2 
followed by 1/ = 1. This described hierarchy of levels 
splitting and sequence of plateaus will be observed inde- 
pendent on cither valleys or spin splitting first. 
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FIG. 3. (Color online) a) Behavior of the resistance at the 
charge neutrality point at fixed Bn and increasing Btot- From 
top to bottom angle and total field: 27° (6 T), 45° (8 T), 63° 
(f2 T), 27° (6 T), 72° (f6 T), 8f° (30 T) b)Behavior of the 
resistance at the charge neutrality point when B has only in 
plane field component(6' = 90°). c) Suggested scheme of the 
spontaneously split valley followed by spin splitting induced 
by B. 



IV. RESISTANCE AT THE CNP IN TILTED 
MAGNETIC FIELD 



In order to clarify the nature of the gapped ground 
state of bilayer graphene and its evolution in magnetic 
field we perform a tilted magnetic field experiment. In 
tilted experiments the total magnetic field (Btot) can be 
separated from its normal (to the plane of the sample) 
component: _B„ = BtotCosB, where 9 is an angle between 
these two vectors (Fig. [3}:). This procedure allows us to 
distinguish between the orbital effect (QHE) and pure 
Zeeman energy, which has to scale with Btot value.— i^i^ 

All measurements presented below were performed at 
a temperature of 1.3 K. The application of the magnetic 
field perpendicular to the sample plane leads to an in- 
crease in the resistance at the CNP, as it is expected for 
a QH transport in the case of broken symmetry states. 
To distinguish between normal component and total B 
we perform a series of experiments with keeping _B„ fixed 
and gradually increasing Btot- As an example, in Fig. [3^) 
we show a change in Rmax at i?„ = 5 T and Btot increas- 
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FIG. 4. (Color online) (a) Change in the R 
peak with total magnetic field Btot- (b) ln(_R, 
tion of Btot at different Bn. The values of B„ from left to 
right 1, 4, 9, 13.7, 17, 21.2, 25 T. (c) The slope of the linear 
fit from Fig. 4b) as a function of normal component B„. 



ing from 6 T up to 30 T for different angles 9. The ac- 
tual maximum of the resistance consists of three peaks: 
highly resistive in the middle {Vg = 1.2 V) and two side 
peaks at the gate voltage at -0.5 and 3 V. The total mag- 
netic field causes decrease in the resistance and the mid- 
dle peak starts splitting into two peaks (or developing 
minimum in resistance at the CNP) when Btot > 6 T for 
studied values of Bn- We observe exactly the same be- 
havior in the experiment when B„ = and applied field 
is in parallel to the graphene membrane: maximum of 
the resistance goes down and develops a local minimum 
at the CNP (Fig. [3J5) . We attribute this change with an 
increase of the total magnetic field. That fact that resis- 
tance changes with Btot indicates that observed effect is 
not a simple quantum localization due to inhomogeneity 
in the sample. 

All three maxima around the CNP decrease in their re- 
sistance in applied parallel B- However, only the middle 
maximum aXVg = 1.2 V shows clear scaling with the total 
magnetic field (Btot) at different tilted angles 9 (Fig.|4ji). 
As one can see in the case of Btot = Bn [9 = 0, black 
curve in Fig. 3^) the resistance keeps on increasing up 
to around 14 T; further increase in magnetic field brings 
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Rmax to lower values (Fig. |4^). Once the non zero angle 
is introduced the common trend for Rmax is a decrease. 

We suggest that the behavior of the middle peak is 
caused by the many-body effect and can be explained by 
the Zeeman splitting closing the spontaneous gap. The 
hierarchy of energy levels is depicted in Fig. [2j;). Once 
B value is high enough the LLL is split in to 4 levels, 
each two-fold degenerate. If we assume that the latter 
degeneracy is spin, then after the appearance of plateau 
associated with filling factor 1^=1 we expect the value 
of the ground state gap A to be lowered by spin splitting 
coupled to Btot- Here we would like to emphasize, that 
we do observe appearance of v = 1 and minimum of 
resistance at the CNP in similar magnetic field Btot > 
7 T. In a simplified way we describe resistance value at 
the CNP point as 

InRmax cx A/(fcT) - g*HBBtot/{kT), (3) 

where g* is an effective g-factor including exchange elec- 
tron interaction and a Landau level broadening i^i^^i'^^ 
The change in ln(i?) versus Btot at fixed i?„ values is 
shown in Fig. . This dependence can be the best de- 
scribed as linear. The slope and y-intercept of the linear 
fit of Fig. |4}d) give the values of A and g*iJ.B- Surpris- 
ingly, these both contributions scale with J5„ component. 
In Fig. I3J;) we show g*fJ,B values versus i3„. Despite the 
fact that the scaling seems like linear, plotting the slope 
as a function of ^/ Bn does seem like fitting as well (not 
shown). A value increases with Bn from 1.4 mcV at 
B„ = 1 T up to 1.7 meV at B„ = 25 T (not shown). This 
A is of the same order as the measured transport gap 
(which can overestimate a real energy gap) and also cor- 
responds to the theoretically predicted gap of 1.5-30 meV 
for the excitonic instabilityJ^ii^i2^ 

In summary, tilted magnetic field experiments show 
that the resistance at the CNP of studied gapped bi- 
layer graphene decreases linearly with the total magnetic 
field component. This points to a many-body effect and 
weak reduction of the gap in applied magnetic field. The 
developed minimum in the resistivity in Fig. [3] can be 
explained by the overlapping of spin-up and spin-down 
levels from the adjacent Landau levels due to Zeeman 
splitting in applied BM- However, from our experiments 
the estimated g* < 0.2, which is very low for spin split- 
ting. In addition, although the resistance decreases in 
parallel field, the Rmax value does not change an order 
of magnitude. This behavior in B is consistent with the 
layer antiferromagnetic state as a ground state of studied 
bilayer samplci^ Since in this state the top and bottom 
layers host spins with opposite orientations, their inter- 
action with applied B can not be described as a simple 
Zeeman splitting. Next to it, our results also open an 
additional question: What is the role of exchange energy 
and level broadening F in LAF state? Naively, scaling 
of g* with Bn can be understood from their depen- 
dence on level broadening F. The F value scales with 
\/Bn, meaning the bigger i?„ the smaller Btot is needed 
to observe level's overlapping. In reality the situation can 



be much more complicated including possible exchange 
mechanisms we do not understand yet. This is also sup- 
ported by the fact that the ground state gap A depends 
on Bn as well. 

Based on these results we suggest a possible scenario 
of symmetry breaking in high quality bilayer graphene 
(Fig.[5J; and Fig. [3};). First splitting is caused by valleys, 
which results in the observed transport gap. An appli- 
cation of magnetic field induces spin splitting of both K 
and K' levels. When B is high enough then the energy 
of spin-up level from K will start approaching the spin- 
down level from K'. The overlapping of the levels will 
cause a decrease in the resistance at the charge neutrality 
point. Since we do observe transport gap in our sample, 
we exclude nematic phase transition. In addition to this, 
the response of the sample in tilted B fits to the LAF 
state. The cause of the valley splitting can be a combi- 
nation of two effects: electron-electron interaction (which 
determines the B field behavior of the middle resistance 
maximum ) and a contamination of the sample surface 
with charged impurities which break inversion symmetry 
(via introduction of electrical field) ^ 



V. CONCLUSIONS 

We report a transport gap of 3 meV in suspended bi- 
layer graphene at 4.2 K, which increases with decreasing 
of temperature. The sequence of appearance of the QHE 
plateaus at the filling factor 1/^2 followed by = 1 
supports a suggestion that the observed gap caused by 
the symmetry breaking. Measurements in the tilted mag- 
netic field indicates that the resistance at the CNP shows 
weak linear decrease with the total magnetic field compo- 
nent. Based on this we propose a possible scenario of the 
symmetry breaking in the investigated bilayer graphene: 
Spontaneous valley splitting at zero magnetic filed fol- 
lowed by the splitting of the spins originating from dif- 
ferent valleys with increasing of B. The gap value and 
weak response of the sample to applied magnetic field 
corresponds to the predicted spin polarized layer antifer- 
romagnetic state as a ground state of the investigated 
sample. The observed non-trivial dependence of the gap 
value from the normal component of B suggests possible 
exchange mechanisms in the system. 
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FIG. 5. (Color online) Electrical scheme of the setup we use 
to perform our measurements. Inset: Scanning electron mi- 
crograph of a typical suspended bilayer membrane in between 
two contacts. 



VI. APPENDIX 



In order to minimize self-heating in graphene at the 
high resistive CNP we used the following scheme (Fig.[5|). 
An AC source maintained a fixed voltage amplitude of 
45 mV (1.87 Hz frequency) across the sample in series 
with 45 Mf2 resistor. The current through the sample 
is monitored by the lock-inl, whose output Ui is propor- 
tional to the current flowing in the circuit (f/i = Ix Ikft). 
Simultaneously, the four probe voltage across the sam- 
ple (C/2) is phase detected by another lock-in2 connected 
through the preamplifier having an input resistance up 
to 100 Then the resistance of the sample is deter- 

mined by i? = Ikfi X U2/U1. The power dissipating in 
the sample is P = C^l/^- Therefore, assuming that max- 
imum V2 is already reached (cx 10 mV), with increasing 
of R the dissipation in the sample will be decreasing. 
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